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We investigate the geometric phases and the Bargmann invariants associated with a multi-level 
quantum systems. In particular, we show that a full set of 'gauge-invariant' objects for an n- 
level system consists of n geometric phases and |(n — l)(n — 2) algebraically independent 4- vertex 
Bargmann invariants. In the process of establishing this result we develop a canonical form for 
U(n) matrices which is useful in its own right. We show that the recently discovered 'off-diagonal' 
geometric phases [ N. Manini and F. Pistolesi, Phys. Rev. Lett. 8, 3067 (2000) ] can be completely 
analysed in terms of the basic building blocks developed in this work. This result liberates the 
^. ' off-diagonal phases from the assumption of adiabaticity used in arriving at them. 
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The notion of geometric phase, though denned originally in the context of adiabatic, unitary and cyclic evolution 
m , has now come to be recognised as a direct consequence of the geometry of the complex Hilbert space and that 
of the associated ray space m . The quantum kinematic [2] picture which has thus emerged provides a much wider 
setting for the notion of the geometric phase by rendering superfluous the various assumptions that attended its 
original discovery (j, and subsequent development |3|-||]. In particular, the requirement of cyclic evolution is no 
longer necessary and it becomes possible to ascribe a geometric phase to any open curve in the Hilbert space which 
has non- orthogonal unit vectors as its end points. Further, following the quantum kinematic approach [pi one is led, 
in a natural way, to the intimate relationship that exists between the geometric phase and the n-vertex Bargmann 
invariants J7fl and that between the geometric phase and Hamilton's theory of turns ||. As an application of this 
approach the Gouy phase (the phase jump experienced by a focussed light beam as it crosses the focus), discovered 
over a hundred years ago, has been shown to be a four vertex Bargmann invariant B. 

In the present work, we develop the quantum kinematic approach for the special case of unitary evolution of an 
n-level system. It turns out that, in the present context, it becomes necessary to introduce Bargmann invariants 
constructed out of two sets of orthonormal basis vectors. We investigate, in detail, their properties and identify and 
construct a full set of gauge- invariant building blocks for the n-level system. We also develop a canonical representation 
of U(n) matrices and bring out its relation to the Bargmann invariants. This representation has recently been shown 
to be extremely useful in parametrizing the Cabibbo-Kobayashi-Maskawa matrix which arises in the context of CP- 
violation in particle physics pQ |. 

As noted above, the geometric phase becomes undefined for those open curves in the Hilbert space which have or- 
thogonal vectors at their ends. A recent work by Manini and Pistolesi [jllj addresses itself precisely to such exceptional 
cases. Employing the original Berry setting of adiabatic evolution for an n-level quantum system, they introduce the 
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concept of off-diagonal geometric phases which can be meaningfully defined in such cases, and showed that these 
off-diagonal geometric phases play an essential role in the interpretation of the findings of a recent experiment ]12fl . 

We study the off-diagonal phases of Manini and Pistolesi within the framework of the quantum kinematic approach 
and show that, in actual fact, this approach, suitably augmented, is robust enough to accommodate the off-diagonal 
phases as well. A brief outline of this work is as follows. In section II we quickly recapitulate the basic ingredients of 
the quantum kinematic approach to geometric phases in a general setting. In section III we specialise the discussion 
to rt-level systems and construct the gauge-invariant objects for this case. These turn out to be n geometric phases 
and a collection of 4- vertex Bargmann invariants. Next we consider the problem of counting, leading to identification 
of a set of independent Bargmann invariants. This necessitates a detailed analysis of the structure of the U(n) 
matrix group. In section IV, we develop a canonical representation for U(n) matrices which is then used, in section 
V, for isolating the independent 4-vertex Bargmann invariants. In section IV, we apply the machinery developed 
in the previous sections to the off-diagonal geometric phases and show how they can be expressed in terms of the 
ordinary geometric phases and the 4-vertex Bargmann invariants, thus liberating these phases from the assumptions 
of adibaticity. Section VII contains our conclusions. 

II. GEOMETRIC PHASES AND BARGMANN INVARIANTS 

We review very briefly the background and ingredients that go into the definition of the quantum geometric phase 
from the kinematic viewpoint, and then the Bargmann invariants and their properties in the generic case. Let Ti be 
the Hilbert space of states of some quantum system, and let B be the set of unit vectors in TL: 

B={yjEH\ U\\ 2 =(yj,yj) = l}cH. (2.1) 

The corresponding ray space (consisting of equivalence classes of unit vectors which differ from one another by phases) 
and projection map are written as 1Z and 7r respectively: 

ir : B — » 7Z = space of unit rays. (2-2) 

To arrive at the concept of geometric phase we begin with parametrised smooth (for our purposes continuous and 
once piecewise differentiable) curves C of unit vectors, which may be pictured as strings lying in B: 

C = {ip(s) eB\ Sl <s<s 2 }cB. (2.3) 

A gauge transformation is a smooth change of phase in a parameter dependent manner at each point of such a curve 
C to lead to another C: 

C = |v'(s) = e ia( - s) tl>(s) | ip(s) G C, ai < s < s 2 } C B. (2.4) 

Then C and C share a common parametrised image curve C in ray space: 

ir[C] = tt[C] = C CTZ. (2.5) 

In general we permit C and even C to be open curves. 

The total, dynamical and geometric phases are then defined as follows as functionals of appropriate arguments: 

PtotPl =arg(V'(si),i/'(s2)), 
PdynPl = Im / ds ( ^( s )' —fo~ ) > 

si 

^ 9 [C}^not[C]-^ dyn [C}. (2.6) 

While the first two phases are functionals of C and do change under a gauge transformation, the geometric phase 
<Pg[0] is gauge-invariant, which explains why it is written as a functional of the ray space curve C. All three phases 
are, however, individually reparametrisation invariant. 

Now we turn to the Bargmann invariants and their 3connection to geometric phases. Given any sequence of n 
vectors ipi , ip2 , • • ■ , ipn in B, the corresponding n- vertex Bargmann invariant is 

A„(^l,V'2)---)V ; n) = (V'l)V'2)(V'2,^3) ■■■ (^n-l,tpn)(^n,tpl)- (2.7) 



Here we assume in the generic case that no two successive vectors in the sequence are mutually orthogonal. It is clear 
that this expression is invariant under cyclic permutations of the t/Vs, and also under independent phase changes of 
the individual vectors. Therefore it is actually a quantity defined at the ray space level. It turns out that the phase 
of A n (tf>i, ip2, ... , ipn) is the geometric phase for suitably constructed closed ray space curves obtained by joining each 
tpj to the next ipj+i (and finally ip n to ip\) by any so-called 'null-phase curve' p|. A null-phase curve is a continuous 
ray-space curve such that for any finite connected portion of it the geometric phase vanishes. That is 'being in phase' 
in the Pancharatnam sense |L3] becomes an equivalence relation on such curves. Examples of null phase curves are 
ray space geodesies (with respect to the well-known Fubini-Study metric |14|Jl5|| ), but the former are a much larger set 
than the latter o. It should be emphasized that whereas the Bargmann invariant (2.7) is defined once its 'vertices', 
namely the projections Tr(ipi),n(ip2), ■ ■ ■ ^(^n) m 7Z-, are given, to interpret its phase as a geometric phase requires 
that we join each ipj to the next ipj+i m some definite manner, namely by some null phase curve, resulting in an 
'n-sided' closed figure in 1Z. 

It can now be seen that as far as phases are concerned, an n-vertex Bargmann invariant for n > 4 can be reduced to 
a product of A3 factors in the generic case [g) , so we can regard the three- vertex Bargmann invariants as the primitive 
ones. For example we have 



A 4 (^i, fa, V>3, ^4) = A 3 (^i, fo, V>3)A 3 (V>i, Vs, ^)/\(i>i, ^s)| 2 , (2.8) 



and more generally 



A n (-0l, ^2, ...,1p n ) = A 3 (V>1, V'2,V'3)A n _i(^i, 1p 3 , 1p4, ■ ■ ■>«)/ 1(^1, ^3)| ■ (2.9) 

Thus the geometric phases of ray space 'triangles', each of whose sides is a null phase curve, are primitive or irreducible 
phases, and all others can be built up from them additively. The purpose in mentioning this is that in the particular 
situation we shall be dealing with later the primitive Bargmann invariants will turn out to be A4's rather than Aa's, 
so that situation will not be generic in the present sense. 

III. GAUGE- INVARIANT PHASES FOR iV-LEVEL SYSTEMS 

We now turn to a study of phases associated with n-level quantum systems. Thus we have an n-dimensional complex 
Hilbert space 7i n describing the pure states of the system. The unit sphere in TL n , and the corresponding space of 
unit rays, will be denoted by £>„ and lZ n respectively. 

If we imagine that a time-dependent Hamiltonian (n x n hermitian matrix) is given, then at each time its complete 
orthonormal set of eigenvectors defines an orthonormal basis for TL n . Assuming there are no degeneracies or level 
crossings, the eigenvalues can be arranged in, say, increasing order; and at each time this basis for TL n is defined upto 
the freedom of phase changes in each basis vector. As time progresses this basis experiences a continuous unitary 
rotation. 

In keeping with the approach of the previous Section, however, we will adopt a kinematic approach here as well 
and not assume any particular Hamiltonian to be given. Thus we imagine that for each value of a parameter s in the 
range s\ < s < S2 we have an orthonormal basis ipj(s),j — 1,2, . ..,n for 7i n ; and as s evolves this basis experiences 
a continuous unitary evolution. Thus we have 

(ipj(s),if! k (s)) = 5 jk , j,k = l,2,...,n; 

^Vj(sMj(s) + =Id, s 1 <s<s 2 . (3.1) 

i=i 

For ease in writing, we shall denote these vectors at the end points s\ and S2 as follows: 

t^.(si) = tpj, ipj(s2) = <f>j. (3.2) 

(The orthonormal vectors ipj here are not to be confused with the arguments of A„ in eqn.(2.7)). Our aim is to obtain 
gauge-invariant expressions and phases in this context. We expect to be able to construct both geometric phases 
Vg[C] for various C, and Bargmann invariants. 

For each value of the index j, as s varies from si to S2, the vector tpji s ) traces out a particular continuous 
parametrised curve Cj in B n : 

C 3 = H'j( s ) £ Bn I 81 < s < s 2 } C B n , j = 1,2, . . . , n. (3.3) 



This curve runs from ipj to (pj. Its image is n[Cj] = Cj C !Z n , and we have n distinct geometric phases: 

Vg [ C j] = notiCj] ~ VdynPj'l' 
"PtotPj] =azg(tl>j,<t>j), 

S 2 

^dynlC,] = Im| ds ^i(«), ^f^J , j = 1, 2, . . . ,n. (3.4) 

«i 

Each of these geometric phases is unchanged under arbitrary alterations in the phase of each ipj (s) at each parameter 
value s. 

We turn next to the construction of Bargmann invariants, the vertices of which are taken from the n initial 
orthonormal vectors ipi, . . . ,ip n and the n final ones <p\ , . . . , <p n . Here we encounter an interesting difference compared 
to the generic situation discussed in the previous Section. Since any two distinct ipj's (and similarly any two distinct 
0j's) are orthogonal, in any Bargmann invariant A„(. . .) an argument ipj must be followed by an argument <p k , which 
must be followed by some ipi, and so on. Similarly if the first argument is some ip, the last one must be some (p. Thus 
in the present context only even order Bargmann invariants A 2 £ survive, a general one being 

A 2 £ (lp n , <j>ki , 1pj 2 , 4>k 2 ■ ■ ■ , 1pj t ,<Pk e ) = (iph , <j>ki ) {4>ki , 1P32 ) {^32 > 4>k 2 ) • • • (ipji i 4>k t ) ((Pkt , V'ji ) (3-5) 

We may now regard the generic case as obtaining when every inner product (ipj, (pk) is nonzero. In this situation we 
find that the primitive Bargmann invariants are A^s rather than Aa's; for instance, 

A6(^j 1 ,<Afc 1 ,V'j2>fe>' ! / ; J3>4 :> fc3) = A 4 (^ , fcl , Vj 2 , <£fc 2 ) A 4 (Vj^ , (p k2 ,ip ja ,(p k3 ) /| {ip ni (pk 2 ) | , (3.6) 

and similarly for higher order A2/S. 

From this discussion it emerges that in the present context the basic gauge-invariant expressions are the n geometric 
phases ip g [Cj] and the various 4- vertex Bargmann invariants A4 : (ipj 1 ,(pk 1 ,ipj 2 ,<pk 2 )- (We are using here quantities 
referring to the entire parameter range s\ < s < $2 and to its end points, and not to any subranges). An important 
problem that now remains is to select out of all possible A^s a maximal set of independent ones as far as phases are 
concerned. For this, as a first step, we turn to an interesting analysis of the structure of the unitary matrix groups 
U(n). 

IV. A CANONICAL REPRESENTATION FOR U(N) MATRICES, COUNTING OF INVARIANT PHASES 

Referring to eqns.(3.1,2) we have a parameter dependent n x n unitary matrix describing the transition from the 
initial orthonormal basis {ipj} for 7i n at s = Si to the moving basis {ipj{s)} at a general s: 

Ms) = (ajk(s)) e U(n), 
djk(s) = (ipj,ipk(s)), Sl < s < s 2 , 
ajk{si) = Sjk- (4.1) 



At s — S2 we write A(s2) = A: 



A = (ojfe), 
djk = (ipj,<f>k)- (4-2) 



The four-vertex Bargmann invariants of the type appearing in eqn.(3.6) are expressions involving products of matrix 
elements of A and their complex conjugates: 

A 4 (lpj,(p k ,1p e ,(p m ) = (ipj , (p k )((p k ,1pi)(lp£, <j> m ){(pm,4>j) 

= a jk a} k a em a* m . (4.3) 

Our problem is to determine how many algebraically independent A4's there are in the generic case in so far as their 
phases are concerned, and to find a convenient enumeration of them. This turns out to be a somewhat intricate 
problem. After some preparation in this Section, the solution will be developed in the next one. 

In working with n x n unitary matrices it is convenient to keep in mind the standard basis in TL n . Then U(n) is 
the group of unitary transformations acting on all n dimensions. For m = 1, 2, . . . , n — 1, we will denote by U(m) the 



unitary group acting on the first m dimensions in 7i„, leaving dimensions m + 1, m + 2, . . . , n unaffected. Then we 
have the inclusion relations (the canonical subgroup chain) 

(7(1) C 17(2) c U(3) ...cU(n-l)C U(n). (4.4) 

General matrices of U(n), U(n — 1), . . . will be written as A n , A n -\, .... In a matrix A m G U(m), the last (n — m) rows 
and columns are trivial, with ones along the diagonal and zeroes elsewhere. (When no confusion is likely to arise, A m 
will also denote an unbordered m x m unitary matrix). 

We will now show by a recursive argument that (almost all) elements A n G U(n) can be expressed uniquely as 
n-fold products 

A n = A„(C)A n _ 1 ( !Z )A n _ 2 (0 . . . A 3 (l)A 2 (a)A 1 ( x ), (4.5) 

where -A n (C) is a special U(n) element determined by an n-component complex unit vector £ G B n ;A n -i(rj) is a 
special U(n ~ 1) element determined by an (n — l)-component complex unit vector r\ G £>n-i; and so on down to 
A^(a) which is a special U(2) element determined by a 2-component complex unit vector a G Bi\ and -Ai(x) is a 
phase factor belonging to U(l). We are led to expect such a representation for A n by the following argument. Any 

vector ( G B n can be carried by a suitable U(n) element into the n vector of the standard basis, (0, 0, ... , 0, 1) T ; 
and the stability group of this vector is the subgroup U(n — 1) G U(n) acting on the first (n — 1) dimensions in 
Ti n . Thus U(n) acts transitively on B n , and this is just the coset space U(n)/U(n — 1). Each coset is thus uniquely 
labelled by some ( G B n . We therefore expect that a general A n G U(n) is expressible as the product A n (()A n -i 
where ( is the last column in A n and A n (Q is a suitably chosen coset representative. Repeating this argument (n— 1) 
times we are led to expect the representation (4.5). The counting of parameters is also just right. Remembering 
that a, (3, . . . , £, r/, £ are complex unit vectors of dimensions 2, 3, . . . , n — 2, n — 1, n and adding the U(l) phase x, the 
number of real independent parameters adds up to n 2 , the dimension of U(n). 

We now present the argument leading to (4.5), yielding in the process the determination of A n {Q). Let a generic 

matrix A n — (ajk) G U(n) be given and let its last (n tn ) column be (,: 



Q, j = l,2,...,n. (4.6) 



Multiplying A n by an A„_i on the right leaves this column unchanged. We choose A„„i so as to bring the n tn row 
of A n to a particularly simple form (for ease in writing we keep using A n and a,jk for the U(n) element obtained at 
each successive stage of the argument): 

a nk = 0, k = 1,2, ...,n-2; 
o-n, n-i — real positive 

= (1-|C«D 1/2 - (4-7) 

The A n _i used here is arbitrary upto an A„_2 factor on its right. Having simplified the n row of A n in this way, 
we can determine all the other elements in the (n — 1) column by imposing orthogonality of rows 1, 2, . . . , n — 1 to 
row n: 

a n ,„-iaj <n -i = -Q 0, j = 1, 2, . . . ,n - 1. (4.8) 

At this point the last two columns and the last row of A n are known in terms of £. 

We next use the freedom in choice of A n -i mentioned above and multiply A n on the right by a suitable A n _2 
(unique upto an A„_3 on its right) to bring the (n — 1) row of A n to a particularly simple form: 

a n -i.k = 0, k = l,2,...,n-3; 

a«-i,n-2 = real positive (4-9) 

Normalising this row gives a n -i,n-2'- 

On.n-1 On-l.n-2 = (l - |C«| 2 - ICn-l| 2 ) V2 • (4.10) 

Next we determine all the remaining elements in the (n — 2) column of A n by imposing orthogonality of rows 
1, 2, . . . ,n — 2 to row (n — 1): 



-1 On-l,n-2 a j,n-2 — - C-lO) J — 1)2, 



(4.11) 



At this point the last three columns and last two rows of A n are known in terms of £. 

This argument can be repeated all the way until we obtain a matrix A n (Q E U(n), all of whose elements are 

determined by the n column £, namely it serves as a coset representative in the coset space U(n)/U(n — 1). (In 
particular the last U{1) element A\{\) is used to make a 2 i real positive). After some algebra we obtain the result 
that the matrix A n (£) = (ftjfc(C)) is uniquely determined by the conditions: 



djj-i(Q = real positive, j = 2,3, 

a jn{() = Cn, j = 1,2, ...,n. 



(4.12) 



Thus A n (() has vanishing matrix elements in the lower left hand triangular portion upto two steps below the main 
diagonal; nonzero matrix elements appear only one step below the main diagonal, and beyond. The explicit expressions 
for the nonzero matrix elements are: 

a j-j-i(0 = Pj-i/pj, 3 = 2,3, . . . ,n; 

a j,k{Q = -Ck+iCj/Pk Pk+i, j < k < n - 1; 
a j«(C) = Cn, 3 = 1,2, ...,n; 

|2 1 ia |2 1 , !/■ ^N 1 / 2 



P, = (ICil 2 + IC 2 | 2 + ... + IC,l 2 ) J 



(1 - IC 



i+il 



IC 



J+2| 



IC 



|2\V2 



Since the quantities /9j obey 



Pi = |Cl| < Pi < P3 < ■ ■ ■ < Pn-l < Pn = 1, 



(4.13) 



(4.14) 



it is evident that this determination of ^4 n (C) goes through with no problems as long as £1 is nonzero, ie., p\ > 0. 

It may be helpful to give the expressions for A 2 (a) € U{2),Az(f3) £ U(3) determined in this way, so as to see the 
general pattern: 



A 2 (a) 



A 3 ((3) 



—ol 2 u\ 

Kl 
M 

P1P2 
Pi I Pi 





Oil 



a 2 



-PIP 



3 PI 



P2 

-/?3*/32 



P2 
P2 



ft 



"2 



(4.15a) 



|/?i| 2 + |ft| 2 + |ft| 2 = l, 



Pi = N> P2 = (1 - |/3 3 | 2 ) 



2\V2 



(4.15b) 



We notice in passing that these are not elements of 5/7(2) and SU(3) respectively. 

Going back to the proof of eqn.(4.5), we see that it can be recursively established; and \, a, (3, . . . , £, 77, C sup- 
ply us with exactly n 2 real independent parameters for A n . Of these, the 2 n ( n ~ 1) independent quantities 
l a i|j \Pi\i \(h\, ■ • • , ICil, IC2I) • • • , Crt — 1 arc of the modulus type, and then there are \n(n + 1) independent phases. 
We can display a general element A n G U(n), (in particular A of eqns.(4.2)) in the selfevident forms 



An = An(C, ■£,£,..., p, a, x) 
= A n (()A n _ 1 , 
A n _! = A„_i(t7,£, ...,P,a,x)- 



(4.16) 



We are now interested in the following operation: suppose we premultiply and post multiply A n by two independent 
diagonal elements of U(n) (a 'gauge transformation' of A n ): 



A n — » A n — D n (0i,6 2 , ■ ■ ■ ,O n ) A n D n (6 1 ,0 2 , . . . ,0 n ) , 
D n (0 1 ,e 2 , ...,*„) = diag (e l6 \e^, . . .,e w ") . (4.17) 

We would like to know: how many independent invariants can we construct out of A n under these transformations, 
how many of them are phases, and how can they be captured through four-vertex Bargmann invariants ? In the 
case of the matrix A of eqn.(4.2) the transformation (4.17) amounts to changing the phase of each ipj and each 4>k 
independently: 



ti 


= e- i0 ^j, 


ftk 


= e<<t>U 


a' jk 


= e<W*) ajk , 



(4.18) 

and we seek an independent set of invariant expressions of the form (4.3). 

First we count the expected numbers of invariants of each kind. The real dimension of U(n) is n 2 . The number of 
independent #'s and 0"s in (4.17) is (2n — 1), because an overall constant phase can be attributed to either the left or 
the right diagonal factor. Therefore the number of real invariants in (n — l) 2 . In the description (4.16) of a general 
A n G U(n), it is clear that under (4.17) every component of each of a, (3, . . . , £ just undergoes a phase change, so the 
quantities |ai|, \Pi\, | /?2 1 , . . . , |Ci|, IC2L ■ ■ ■ j ICn-il are \n(n — 1) real independent modulus type invariants. Thus there 
must be a balance oi \{n~ l)(n — 2) real independent phase type invariants. This is in agreement with known results 

We shall describe in the next Section a recursive procedure by which we can pick out ^(n — l)(n — 2) algebraically 
independent four- vector Bargmann invariants whose phases are the expected phase invariants associated with a general 
U{n) matrix. 

V. DETERMINATION OF INDEPENDENT BARGMANN INVARIANTS 

We describe first how, in a recursive manner, we can isolate the expected ^(n — l)(n — 2) independent gauge- invariant 
phases for a generic A n £ U{n) using the parametrisation (4.5), and then turn to the choice of an equal number of 
independent primitive Bargmann invariants A4. 

We begin with eqns.(4.5,4.16), 

A„_i = A„_i(t7)A„_ 2 (0 . • • A 2 (a)A 1 ( X ) 

= A„_i(?7,£, ...,a,x), (5.1) 

apply diagonal matrices on the left and on the right as in eqn.(4.17), and trace the changes that occur in £ and in 
A n - X : 

A' n = D n (6i,6 2 , ■ ■ ■ ,0 n )A n D n {9'i,0' 2 , .. . ,9' n ) 

= A n {OA' n _ 1 . (5.2) 

Our aim is to compute £ and A' n _ 1 . Since the D n factors are quite elementary this can be carried through as follows: 

A' n = £>„(#!, 2 , . . .,e n )A n {QA n - X D n (0' l7 e' 2 , ...,6' n ) 

= D n {e Xi e 2 , ..., 9 n )An(gAn-i D n (o, o, . . . , o, e' n ) D n (e[,e 2 , . . . , e' n _ x , o) 

= D n (e!,9 2 , . . . , e n )A n (QD n (0, . . . , 0, 6' n ) A n _ t D n _ x fax, . . . , C-i) ■ (5-3) 

The product of the first three factors simplifies: 

D n (6 1 ,6 2 ,...,6 n ) A n (gD n (0,...,0,6' n ) 

= D n (0 1 ,e 2 , . . .,e n )(a jk (0)D n (0,0, ...,e' n ) 
= (MC)), 

b jk (0 = e^o ifc (C), k = 1, 2, . . . ,n - 1. (5.4) 



Here the matrix elements CLjk(Q are given in eqn.(4.13), and for simplicity the 9 and 9' dependences of bjk are left 
implicit. In particular, as in eqn.(4.12), 

MC) = o, k 

while 



l,2,...,j-2; j=3,4,...,n, 



(5.5) 



bjj-i(C) = e»>a j , j - 1 (g 



j Pj-i/pj, j = 2,3,. 



(5.6) 



Thus the matrix (&jfc(C)) would have been (ctjfc(C')) except for the fact that the elements 6j,j-i(C) J us t below the 
main diagonal are not real positive but carry phases. But this can be easily taken care of by extracting a suitably 
chosen diagonal matrix on the right: 



(b jk (Q) = (a 3k (C)) D n (6 2 ,9 3 , . . . ,8 n ,0). 



(5.7) 



The point is that, according to the statement accompanying eqns.(4.12), after removal of this diagonal factor what 
remains is necessarily A n ((') = (ajfc(C'))- Combining the above steps we get: 



A n — j4 k (C )A»-i 

= (b jk (())A n _ 1 D n _ 1 (d' 1 ,e' 2 ,...,e , n _ 1 ) 

= A n (C)D n (e 2 ,e 3 , ...,e n , o)A n _ 1 d„_i (e[,e 2 , . . . , < 

so the changes induced in £ and in A n _\ by the gauge transformation (5.2) are: 



c; = ^ +e -)o, i = 1,2, 



, n; 



A n _\ — D n -i(9 2 



9 n )A n -i D n _i \9 l ,9 2 . 



(5.8) 



(5.9a) 
(5.9b) 



We see from the structure of this result that we can tackle our problem recursively: The gauge transformation(5.2) at 
the U(n) level translates into the change ( — > £' given by eqn.(5.9a) and a gauge transformation A n _\ —* A' n _ 1 at the 
U(n — 1) level given by eqn.(5.9b). Therefore all gauge-invariant expressions that exist at the A n _i or U(n — 1) level 
survive when we move from U(n — 1) to U(n), and in addition as the vector C, G B n becomes available, new invariant 
phases involving Q can be constructed. The number of the latter can be immediately computed: it is the difference 
between ^(n — l)(n — 2) and |(n — 2)(n — 3), namely the difference between the numbers of gauge-invariant phases 
at the U(n) and the U(n — 1) levels, and this is (n — 2). Therefore the number of new independent phase invariants 
involving A n (Q, ie., (, in an essential way must be (n — 2). These can now be isolated or explicitly constructed as 
follows. 

From eqn.(5.9) we notice that 9\ and 9' n appears only in the transformation law for £, not for A n -\. Therefore we 
first form the (n — 1) independent combinations C*Cj+i to eliminate 9' n completely: 



ac 



j w+i 



-i(6j-6j+i)f* 



CjO+1' 3 



1,2, 



(5.10) 



Here 9\ occurs only in the transformation law for C1C2, being absent as we just mentioned in the law for A n _\. Next 
we notice that the phases 9[, 6' 2 , . . . , 9' n _ 1 , involved in A' n _ 1 , arc completely absent in the transformation law (5.10) 

of QQ+i- Let us therefore look at the (n — 1) column, say, of A n _i, which as is evident from eqn.(5.1) is just the 
(n — 1) component complex unit vector rj G S„_i: 



A n -1 = 



( 







m o\ 

V2 
Vn-1 

1) 



The "earlier" columns of A n _\ are more complicated, as is clear from the structure of A n _\ in eqn.(5.1) 
eqn.(5.9b) we can read off the transformation law for the t/s under the gauge transformation (5.2): 



(5.11) 



From 



v = p^+i+eUO 



Vj, j = 1,2, 



,n — 1. 



(5.12) 



To eliminate 0' n _ 1 we form the (n — 2) combinations fJjVj+i which transform thus: 

VjVj+i ^ e'^'^'^Wj+i' i = l,2,...,n-2. (5.13) 

Comparing cqns.(5.10) and (5.13) we immediately obtain the expected (n — 2) independent (phase-type) invariants 
involving £ £ B n in an essential manner, namely they can be taken to be the complex quantities 

VjVj+i(j+i^+ 2 > 3 = l > 2 ' • • ■ ' n ~ 2 - ( 5 - 14 ) 

By recursion the complete set of 5 {n — l)(w — 2) independent phase type invariants that can be formed from a generic 
matrix A n G U{n) can be written down in terms of the canonical parametrisation (4.5) for A n , and the list reads: 

a J a j+i0j+i0j+2> j = 1; 
/3j/?j+i7j+i7j+2, .7 = 1,2; 



O^+i^j+1%+2. i = 1, 2, • • • ,n - 3; 

W* + iC*+iG+2, j = l,2,...,n-2. (5.15) 

While we have here an explicit solution to our problem, the difficulty is that these invariants are not directly 
expressed in terms of the matrix elements of A n = (ajk) £ U(n). It is true that in our parametrisation £ is the last, 

n , column of A n ; but the previous, (n — 1) , column involves both r\ and £; the (n — 2) column involves £,77 
and C, and so on. The task that remains is to see how to translate the expressions (5.15), as far as their phases are 
concerned, into an algebraically equivalent set of ^(n — l)(n — 2) expressions formed as simply as possible out of the 
matrix elements of A n . We turn to this now, bringing in the 4- vertex Bargmann invariants of A n . 

As indicated in eqn.(4.3), a general 4- vertex Bargmann invariant requires the choice of some two rows, say j and £ 
with j < £, and some two columns, say k and m with k < m, and use of the four matrix elements at their intersections: 

A-jikm = ajkal k a£ m a* m . (5.16) 

But as far as phases go, we can see in a step by step manner that a general Ajikm reduces to a product of factors of 
the simpler form 

Ajk = Ajvy+i^fc+i (5-17) 

involving some two adjacent rows and some two adjacent columns. This is to be understood modulo real positive 
definite factors coming from the squared moduli of some of the matrix elements of A n . The two "recursion formulae" 
that help us achieve this simplification are: 

Aj£km — Aj£km-1 A^ m _i m /|a^ m _ia^ m _i| , 

= Aj£-ikm A£-iekm/\ae-i,ka£-i, m \ ■ (5.18) 

It therefore suffices to work with the (n — l) 2 expressions 

Ajk = ajk a* j+ i,k a j+i,k+i a* hk+1 , j, k = 1,2, . . . ,n - 1, (5.19) 



and their phases. Our goal now is to (at least in principle and in the generic situation) express (the phases of) the 
2(71 — l)( n ~ 2) complex invariants (5.15) in terms of (the phases of) the (n — l) 2 complex invariants (5.19). (In this 
process any number of real positive factors may intervene). We already have here an indication that the (n — l) 2 
expressions (5.19) (more exactly their phases) cannot all be independent, the number of independent ones being only 
i(n — 1)(" — 2). It will turn out, as we indicate below, that these may be taken to be the Ajk for j < k < n — 1. 
Again the proof is recursive in nature. 

Consider the (n — 1) invariants (5.14) that get added to all previous ones when we make the transition U(n — 1) — » 
U(n) and bring in the vector £ £ B n . Instead of being expressed in terms of £ and r\ 6 B n ~\, we now show that they 

can be equally well expressed in terms of C, and the penultimate, ie (n — 1) , column of the complete U(n) matrix 
A n . Let us denote this column vector by w £ B n ] it is orthogonal to £. As noted earlier, it is easily determined in 
terms of £ and 77 or, more conveniently for our purpose, rj is expressible in terms of w and £. Starting with 



A, 




and transposing A n {() we get 



A„(C)A„_i(? ? ) J 4„_2, 



A„_i(?7)A„_2 =A„(C)U r , 



(5.20) 



(5.21) 



Since the factor A n -2 does not affect the last two columns on both sides, we can use the matrix elements (4.12,13) of 
A n (() to obtain: 



% = X! a fej(c)*wfe 



fc=l,2, 

Pj + 1 



Wj+1 



c 



J'+l 



Pjpj+l 



Y^ Q w ^ i = 1 ' 2 ' 



, n — 1. 



(5.22) 



fc=i 



The gauge transformation laws of £ and 77 are given in eqns.(5.9a,12), while that of w is seen from eqn.(5.2) to be 

Wj — >e i ( 6i+6 ' n - 1 ')w j , j = l,2,...,n. (5.23) 

Naturally the relations (5.22) are consistent with these transformation laws. The combinations of 77 and £ needed 
in (5.14) are ?7jC*+i ^ or 3 = 1>2, ...,n — 1. We see from eqns.(5.22) that they are real linear combinations of 

WlCl, W2C2 > ■ • • ' W n-lQ-l-,W n Q 



VjC+i 



10+1,2 £c>, + -^ +1 c + , 



pjpj+i 



fc=i 



pj+i 



(5.24) 



Using both the orthogonality of w and £, and the reduction process (5.18) for A^s formed out of the last two column 
of A n , it is now clear that the set of complex invariants (5.14) can be replaced by the following set of (n — 2) Azi's: 



Aj,„_i = Wj.CjWj^Q+i, j = 1,2, . . . , n - 2. 



(5.25) 



The known algebraic independence of the set (5.14) implies a similar independence of these A^s. 

To tackle the next set of (n — 3) invariants £j?7* +1 £* +1 r7 : ,+2 for j = 1, 2, . . . ,n — 3 in the list (5.15), we must bring 

in the (n — 2) column of the matrix A n . Denote this by v G B n so that 



/ 



A n 



Vl W\ Cl \ 

V2 W 2 C2 

V n W n Cn J 



(5.26) 



Analogous to cqn.(5.21) we now have 



A n _ 2 {Z)A n _ 3 = A n -l{^ A n (C) ] A n , 



(5.27) 



from where we get expressions for £j in terms of £,77 arL d v. This is naturally more complicated than eqn.(5.22) at 
the previous stage. In place of the real positive factors pj defined in terms of ( in eqn.(4.13), we now have similarly 

defined factors <jj in terms of r\ occurring in the elements of A n -\{r\). The result of comparing the (n — 2) columns 
of both sides of eqn.(5.27) is: 
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J'+l 

(TjPj + 1 (T, 






-»?fcVfe+i 



CTj+lPj+2 (Tj+lpj+lpj+2 ^ CTj'^+1 " Pfe+1 

+ ^ ± ^E^ ±1 EC;^ i = l,2,...,n-2. (5.28) 

cr J <r i+ 1 ^ PkPk+i fri 

Here next we can use eqn.(5.22) to go from r\ to w. Then we form the expressions £,jf]* +1 and step by step work our 
way up to the invariants ^jVj+it.j+iVj+2- We can then see that apart from various real factors we encounter A^s 
involving v'a and id's, w's and £'s and id's and £'s. Using the reduction rules (5.18) the v — C combinations can be 
eliminated in favour of the other two types. It is now clear that apart from the A Jjn _i in eqn.(5.25) which appeared 
at the previous stage, the new quantities that come in now are A :/j „_2. But we know in advance that at this stage 
only (n — 3) new independent invariants are available. As all the rows of A„ are on equal footing, we conclude that 
the new A^s to be added now to the previous Aj :n _i may be taken to be 

Aj>_2 = VjWjVj +1 Wj+i, j = 1, 2, . . . , (n - 3). (5.29) 

In this manner one sees recursively that the h (n — 1) (n — 2) independent gauge-invariant phases in a general matrix 
A n £ U(n) are the Bargmann invariants Ajk for j < k < n — 1. In any case such a choice is permitted. However the 
actual algebraic expression of a general Ajk in terms of this special subset may be rather involved, so one may freely 
use all Ajk in constructing interesting gauge-invariant expressions with various properties. 

The upshot of these considerations is that the naturally available gauge-invariant phases for the continuous unitary 
evolution of an n-level quantum system, barring degeneracies and level-crossings, are n geometric phases <p g [Cj] as 
defined in eqn.(3.4), and the (n — l) 2 primitive four-vertex Bargmann invariants Ajk of eqn.(5.19); of the latter, only 
the h (n — 1) (n — 2) Ajk 's for j < k < n — 1 are independent. Any composite expression formed out of these ingredients 
is of course also invariant. 

VI. OFF-DIAGONAL GEOMETRIC PHASES 

It is evident from the definitions (2.6) that while the dynamical phase (pj vn [C] is always numerically well-defined 
once the parametrised curve C is given, the total phase y^otP] ^ s on ly defined modulo 27r, and moreover is undefined 
if the vectors ip( s i) an( i V'( s 2) at the end points of C are orthogonal. These properties naturally carry over to the 
geometric phase y s [C]: only defined modulo 2ir, undefined when </^otP] ^ s undefined. The Bargmann invariants (2.7) 
too share these problems of definition as far as their phases are concerned, which explains the limitation to generic 
situations. 

Recently a very interesting attempt to define so-called off-diagonal geometric phases has been made to cover just 
these exceptional or problematic situations |TJ]. Specifically the idea is to set up gauge- invariant phases associated 
with the unitary evolution of an n- level quantum system, which remain well-defined even when one of the eigenvectors 

of the Hamiltonian at a final time t-z, say the k one, happens to coincide with the j eigenvector at the initial 
time ti, with j ^ k. In this situation, as ipj(ti) and V'fc(fe) are the same upto a phase, both the geometric phases 
ip g [Cj] and ip g [Ck] become undefined since the inner products ("0^ (ti ) , "0j (*2 ) ) and (tpk(ti),ipk(t2)) vanish. We shall 
now briefly recall the basic quantities introduced in this new approach, and then show that the usual geometric phases 
and Bargmann invariants as defined earlier can completely handle the new situation. It is just that they must be 
put together in such combinations so that the potentially undefined factors in each precisely cancel one another in 
exceptional situations. 

The notation for the evolution of an n level quantum system is as given in Section 3. The quantities defined in the 
off-diagonal geometric phases method, when expressed in our notations, are: 

Ij =exp j-i<^ dyn [Cj]j, j = 1,2, ...,n; (6.1a) 

a-jk = cxp |i arg(-0j, 4>k) ~ ^dynI Cfc ]} > J ^ fc ; (6.1b) 

Ijk = <?jk (Tkj, j ¥" fc ; (6-lc) 

7j- =exp{i(f g [Cj]}, j = l,2,...,n. (6. Id) 

Of these, Ij and ajk are not gauge-invariant, but jjk and jj are gauge- invariant. In case for some j ^ k we have 
\{ipj,4>k)\ — 1, it is clear that both <p g [Cj] and ip g [Ck\ become undefined, but the "off-diagonal" quantity ^jk remains 
well- defined. 
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The two-state or two-index quantity 7^ has been generalised to a multi-index quantity of order i as follows: 

r YjlJ2---je ~ <7 jlJ2 <J J2J3 • ■ ■ a 3e-ije' 7 3eji> \°-^) 

and this again is gauge- invariant. 

We can now see that all these newly introduced gauge-invariant off-diagonal quantities ljk,lj 1 j 2 ...j e are actually 
expressible completely in terms of the geometric phases and Bargmann invariants for the n-level system, in carefully 
chosen combinations: 

7j fc = cxp { % arg A 4 (ipj , </> fe , ^ fc , ^ ) +iip g [C 3 } + iip g [C k ] } , 
lhh-ii = c xp{ i & T E A 2e(<i)j 1 ,ipj 1 ,(j)j 2 ,i>j 2 ,---(j) M ,ip M )+itPg{C J1 }+itp g {C : j 2 }+ ... +iip g [C M }}. (6.3) 

In the case of "fjk, for example, we see that when \{i/)j , <j)k)\ = 1 an d (p g [Cj\, (p g [Ck\ become undefined because the 
total phases ^totPj] an d ^totPfe] are undefined, there are compensating factors from A^tpj, (f>k, xpk, 4>j) that precisely 
cancel these parts of the individual geometric phases, so that jjk remains unambiguous. The mechanism is similar in 
the case of the higher order expressions Jj^.-.ji 

It has been shown that r Yj 1 j 2 ...j e f° r ^ > 4 can be reduced to the expressions with 1 = 2 and I = 3, so these are the 
primitive ones. Among these, we can limit ourselves to choices obeying j\ < ji when £ = 2 and j\ = 1 < ji < j'3 when 
(. = 3, in counting independent quantities. However the upshot of our analysis is that we can always work with just 
the geometric phases <p g [Cj] and the independent A^s listed in the previous Section (but for convenience employ all 
the Ajfc if necessary). All gauge- invariant quantities can be built up out of them, so that conceptually the off-diagonal 
geomeric phases are constructed out of previously known familiar building blocks. 

VII. CONCLUDING REMARKS 

We have carried out a complete analysis of the gauge-invariant objects for n-level quantum systems. This entails 
introduction of Bargmann invariants defined over two sets of orthonormal basis vectors, demonstration that the 
primitive Bargmann invariants are 4-vertex Bargmann invariants, and finally the identification of an algebraically 
independent set of 4-vertex Bargmann invariants which turn out to be (n — l)(n — 2)/2 in number. In the process of 
achieving this task we developed a canonical form for U(n) matrices in terms of a sequence of complex unit vectors 
of dimensions n, n — 1, . . . , 1 which may be useful in other contexts as well. Indeed, this form has already found 
application in parametrising the CKM matrices which arise in the context of CP-violation in particle physics. The 
gauge-invariant building blocks constructed here are shown to provide a complete quantum kinematic picture of the 
recently discovered off-diagonal phases. The usefulness of the off-diagonal phases is thus extended far beyond the 
restrictive framework of adiabatic evolution. This reinforces the view that the Bargmann invariants and the traditional 
geometric phases, and suitably constructed combinations of them, suffice in answering all interesting questions in this 
domain. 
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